Cosmic Evolution in Brans-Dicke Chameleon Cosmology 
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Abstract: We have investigated the Brans-Dicke Chameleon theory of gravity and obtained 
exact solutions of the scale factor a(t), scalar field an arbitrary function which interact 

with the matter Lagrangian in the action of the Brans-Dicke Chameleon theory and potential V{(j>) 
for different epochs of the cosmic evolution. We plot the functions a(t), 4>{t), f(t) and V{cj>) for 
different values of the Brans-Dicke parameter. In our models, there is no accelerating solution, only 
decelerating one with q > Q. The physical cosmological distances have been investigated carefully. 
&X)' Further the statefinder parameters pair and deceleration parameter are discussed. 

< 
s, 

Q-i The Brans-Dicke (BD) theory of gravity defined by a scalar field (f) and a constant coupling function w [l| , is perhaps 
(— ( ' the most natural extension of general relativity (GR) which is obtained in the limit of w — > c» and ()>= constant [2M4!| . 
^ , An imperative property of the BD theory of gravity is that it yields simple expanding solutions [5] for scalar field 
bJO and scale factor a{t) which are well-matched with solar system experiments [GHsj. There are many works on 

interesting physical aspects of the BD theory i9|]. 
O hi a recent paper the dynamics of the cosmic evolution has been investigated in the formalism of generalized BD 

theory of gravity (where uj is not a constant but a function of ^1^. There a consistent solution of the generalized 
BD equations of motion based on simple power-law temporal behavior of a, (j) and u) is obtained. 
' In the present paper, a BD theory in which there is a non-minimal coupling between the scalar field and the 
I matter field is considered. Thereby the action and the field equations are modified due to the coupling of the scalar 

field with the matter. In the literature such type of scalar field usually called chameleon field [11]. We note that a 
04 , chameleon field requires an strict form of its potential in order to avoid the appearance of fifth forces or violations 
^ of the Equivalence Principle, so that a chameleon model has to present a realistic potential. This is due to the fact 
that the physical properties of the scalar field, such as its mass, depend delicately on the environment. Moreover, in 
high density regions, the chameleon mix together with its environment and becomes essentially invisible to searches 
' for Equivalence Principle violation and fifth force [ll| . Further more, it was shown that in the presence of chameleon 
field, all existing constraints from planetary orbits, such as those from lunar laser ranging are easily satisfied [lllfl^. 
The explanation is that the chameleon-mediated force between two large objects, such as the Earth and the Sun, is 
much weaker than one would bluntly expect. In particular, it was shown that the deviations from Newtonian gravity 
due to the chameleon field of the Earth are suppressed by nine orders of magnitude by the thin-shell effect [l^] • Some 
other studies on the chameleon gravity have been investigated in [l3l - [T5| . 

Our work differs from that of Ref. [loj in that we assume a non-minimal coupling between the scalar field and 
the matter field, besides, we are taking as a constant and not a function of </'(i). Our aim is to investigate the 
;h ■ dynamics of the potential V{(t)) and the function /(0) (these are defined in the next section) alongside the scale factor 
a{t) and the scalar field (j){t). In other words we are interested to find exact solutions of a{t), (f>{t), f{(f>) and V{(j)) for 
different epochs of the cosmic evolution. The stability analysis of the solutions is important. But it will be beyond 
the scope of our paper. Such analysis can be done using the similar methodology, as it has been done in [lg|, by 
defining dimensionless variables. The stability of BD with Chameleon field has been done already in [13] for the same 
system of equations as we have. 

The statefinder parameters pair {r, s}, allows one to explore the properties of dark energy (DE) independent of the 
model [l^. It has been used to distinguish fiat models of the DE. Recently this pair has been evaluated for different 
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models [l9|-|30|. In the framework of Brans-Dicke Chameleon cosmology, the statefinder parameters are studied in 

Plan of the paper is as follows. In the next section we give the basic equations of cosmic evolution with chameleon 
scalar field. In the section III we discuss our model and obtain exact solution of a{t), (l>{t), f{(f)) and V{(t)) for (A) 
radiation dominated era, (B) dust fluid era and (C) vacuum energy dominated era. In the section IV the statefinder 
parameters and the deceleration parameter are investigated. The cosmological distances have been discussed in section 
V. Finally we conclude our discussion in section VI. 



II. COSMIC EVOLUTION WITH CHAMELEON SCALAR FIELD 



We begin with the BD chameleon theory in which the scalar field is coupled non-minimally to the matter field via 
the action [s^l 

S = j d^X^ (^^R - ^gf^-d^(pd,(^ ~ + , (1) 

where R is the Ricci scalar curvature, is the BD scalar field with a potential V {(()). The chameleon field (j) is 
non-minimally coupled to gravity, lj is the dimensionless BD parameter. The last term in the action indicates the 
interaction between the matter Lagrangian Lm and some arbitrary function /(</>) of the BD scalar field. In the limiting 
case /((/)) = 1, we obtain the standard BD theory. 

The gravitational field equations derived from the action ([iJ with respect to the metric is 

1 o _ /('^) 



+ ^[(t>fj^;,y - 9^L,^0(f>] - 9^l^l^^-^■ (2) 

where T^i, represents the stress-energy tensor for the fluid filling the spacetime which is represented by the perfect 
fiuid 

Tfj^i^ = (p + p)Uf^u„ + pg^^, (3) 

where p and p are the energy density and pressure of the perfect fluid which we assume to be a mixture of different 
kinds of matters. Actually there is only chameleon field, whose pressure and density are given by a perfect fluid 
stress tensor. For different state parameters, this fluid behaves differently, like matter, radiation or DE. Also u^^ is 
the four- vector velocity of the fluid satisfying u^u^ = — 1. The Klein-Gordon equation (or the wave equation) for the 
scalar field is 

where T is the trace of ([3|)and □ = V^V^ in which the operator represents covariant derivative. The homogeneous 
and isotropic Friedmann-Robertson- Walker (FRW) universe is described by the metric 

ds^ = -df + a\t) [j^^ + r^dn"^ , (5) 

where a{t) is the scale factor, and k = —1, 0, -1-1 corresponds to open, flat, and closed universes, respectively. Variation 
of action ([T|) with respect to metric ([5]) for a flat k = universe filled with perfect fluid yields the following field 
equations 

(6) 



where H — d/a is the Hubble parameter. Here, a dot indicates differentiation with respect to the cosmic time t. The 
dynamical equation (energy conservation) for the scalar field is 
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Similarly the energy conservation for the cosmic fluid is 

p + 'iH{p + p)=Q. (9) 
We shall use the equation of state (EoS) for the fluid p = jp, thus ^ yields 

p = Ca-3(i+7). (10) 

III. OUR MODEL 

Observational data of SN la suggests that theoretical models based on power-law forms of the Chameleon potential 
and scalar function f{(/)) are consistent with the data 'si']. Hence we shall follow the procedure of and will 
obtain solution of the above dynamical equations © to by assuming power law dependence on time for a{t), (t>(t), 
/(</>) and Vi^). 

ait)^ao[-^y, (11) 
Ht)^<l>o(rY, (12) 



fm))-m''^h{lT\ (13) 



vm))-m'' = vo{iy^. (14) 



.to' 

Note that w is a constant BD parameter and ag, (j)Q, fo and Vo are also constants. Also notice that dynamical system 
is a closed system i.e. four differential equations (O to (jO)) for four unknown parameters {a, /3, m, n) to be determined. 
a must be positive for an expanding Universe while other parameters are free. Using in pO|) . we get 



.to/ V to 

Using the set of ansatz functions(ll-15) we get 



/t\-3a(l+7) /tx-3Q(l+7) 

P = Po(-j , p = MtJ ■ ^^^^ 



--13^ +al3 ^0, (16) 
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/oPo = -y, (17) 

(n- m)^ = 3a(l + 7). (18) 

We now proceed to check the consistency of above equations with Eq. (8) which is the wave equation for the scalar 
field 4>{t). Using eqs. (11-15) Eq.(8) reduces to 

/3(/3 - 1 + 3a) = 0, (19) 
(2-m)Uo = Po/o(2-n)(l- 37), (20) 
(TO-n)^ = -3a(l-|-7). (21) 

The above Eq. (19) implies that wc can have /3 = or /? = 1 — 3a. The first one results that 

Ht) = 00, (22) 
/(0) = fo, (23) 
V{(t>) = Vo. (24) 

The above set of functions or their equivalent value f5 — in (16) gives us a = 0. But we remove (3 — in our toy 
model we must take a > 0. Thus in the present work we discard this special case and limit ourselves only to 

/3 = 1 - 3a (25) 

Now we analyze different cosmological epochs. That is the cases of 7 = i, 7 = 0. We remove completely the case 
with 7 = — 1, since as we pointed it previously , in our models, there is no accelerating solution, only decelerating one 
with q > 0. 
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FIG. 1: Variation of scale factor °+^*^ for different values of the BD parameter w. 



A. Radiation dominated (7 = 1/3) 
In this case Eq. (20)gives m = 2. Substituting /3 = 1 — 3a in (16) we obtain 



3(6u + 1) ± Vfe7+9 

10 I n.. ' (26) 



12 + 9c^ 

1 T V6a; + 9 



Now from (18) we obtain 



Thus we have 



^ _ 14+12i^±2v'6u>+9 



(27) 



14 + 12a; ± 2V&7+9 
n = , (28) 



, 3(i^ + 1)±v/6^:J+9 

a(t)=ao(-) , (29) 

^(i)=0o(^)'^', (30) 



/(</>(*)) = /o(-) , (31) 



, 2 lq= V6aJ+9 

mt)) = yo(^) . (32) 

Below we plot some figures include the time behavior of the set of functions a{t), (pit), f{(j){t)), V{4>{t)) for some large 
values of the parameter oj. The first two figures show the accelerated expansion of Universe for some large values of 

the BD parameter. 

The next two figure show the two possible values of the scalar field for a sample of large BD parameter. 
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FIG. 2: Variation of scale factor for different values of the BD parameter ui. 
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FIG. 3: Variation of the one possible value of the scalar field for a sample of large BD parameter w. 



Wc obtain that the radiation dominant toy model is a simple square power model. Here we plot the potential 
function as a function of the scalar field for some values of the (po^Vo- 

Remind that for the limit of the very large BD parameter a; which we expect that the theory must has a GR limit 
, the exponents and the functions are 
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n ~ 0. 



(33) 

(34) 
(35) 



aGR{t) = "o(^) 



(36) 
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FIG. 4: Variation of the another possible value of the scalar field for a sample of large BD parameter u). 



f+it) 



25000- 



20000- 



15000- 



10000- 



5000' 





/ 




/ ■■■/ 




/ .■■/ 




/ 




/ 








.■>' 







2 4 6 







8 10 



I (0=10 o>=100 co=1000| 

FIG. 5: Variation of for a sample of large BD parameter w. 



vm)) - Vo. 

The plot of the scale factor for the GR limit is Fig. (10) 
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FIG. 6: Variation of f-{t) for a sample of large BD parameter u. 
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FIG. 7: Variation of V+(t) for a sample of large BD parameter w. 



B. Dust fluid (7 = 0) 

This case from (20,21) we've 



m+|=3 (40) 



9(a; + 1) ± 3v/fe7+9 

(to - n) = /^-^ (41) 

1 =F V 6w + 9 
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FIG. 8: Variation of V- {t) for a sample of large BD parameter w. 
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FIG. 9: Variation of V^(</>) for a sample of the parameter o = . 



from (40,41) we obtain 



n = 



2(3a; + 4) 



m 



Thus we have 



1 T \/6a; + 9 
-(1 + 3cj)=f3V6w + 9 
1 T V6w + 9 



/Wi)) = /o(^'' 



(42) 

(43) 

(44) 
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FIG. 10: Variation of scale factor in GR limit 



Vm)) = V,[-) '^'^ . (45) 



for large values of the BD parameter 7, we have V{4){tj) = Vq 



IV. STATEFINDER PARAMETERS AND DECELERATION PARAMETER 



The statefinder parameters r and s depends on the third and second derivatives of the scale factor a, just as 
the dependence of the Hubble parameter H and the deceleration parameter q on its first and second derivatives 
respectively. 

The deceleration parameter is defined as 

g = (46) 



The statefinder parameters are [I^ 

'aa^ r — 1 



s = 



3(9-^) 

For Eq. (29), the above three parameters take the form 



(47) 



3 + 2a; =F V9 + 6w 
9± = , 



(ip21 =F ISw + V9 + 6w)(=f9 T 6a; + V9 + 6w) 
^ (3 + 3a; + V9T6ZJ)2 ' 



6 + 6cj T \/9 + 6^ 
^± = 3^^ ■ 



V. COSMOLOGICAL DISTANCES 



In this section we discuss the time dependent cosmological distances of the models presented in sections A,B. 
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FIG. 11: Variation of the deceleration parameter q+,-- The line is q+ and dot denotes g_ 
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FIG. 12: Variation of the parameter r+,_. The line is r+ and dash denotes r_ 



A. Lookback time 



If a photon is emitted by a source at the instant t and received at time ioj then the photon travel time or the 
lookback time 1 — 1^ is defined by 



dz 



(49) 



lo H{z){l + zy 

where oq is the present value of the scale factor of the Universe. If a photon emitted by a source and received by 




FIG. 13: Variation of the parameter s+,-. The line is s+ and dash denotes s- 
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FIG. 14: The behavior of the lookback time as a function of the redshift for values u). 



an observer at time to then the proper distance between them is defined by 

t-to = {l + z)-^. 



(50) 



The Figure- 14 shows the behavior of the lookback time as a function of the redshift for value different values of the 
u) and for the plus sign in (26). 
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FIG. 15: The behavior of the proper distance as a function of the redshift for values uj. 



B. Proper distance 

If a photon emitted by a source and received by an observer at time to then the proper distance between them is 
defined by 



= Oo 



da 



dz 



, am Jo H^{z){l + zy 



which for (29) simphfies to 



d = — (1 



2a' {l + z)i 



(51) 



(52) 



C. Luminosity distance 

If L be the total energy emitted by the source per unit time and i be the apparent luminosity of the object then 
the luminosity distance dL = (3^) ^ evolves as 

dL = d{l + z). (53) 
The figure (16) shows the variation of the luminosity distance as a function of the exponent lo. 



D. Angular diameter 

The angular diameter distance dA is simplified to 

dA = dL{l + z)-"^ . 

The Figure- 17 shows the variation of the angular diameter as a function of the exponent < n < 1. 



(54) 
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FIG. 16: The behavior of the luminosity distance as a function of the redshift for values w. 
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FIG. 17: The behavior of the angular diameter as a function of the redshift for values w. 
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VI. CONCLUSION 



To recapitulate, we studied the Brans-Dicke Chameleon cosmology by obtaining exact solutions of the scale factor 
a{t), scalar field (/'(t), the potential V{(j)) and the arbitrary function /((/>) for different epochs of the cosmic evolution. 
This analysis was performed by using ansatz for the above parameters. These were motivated since earlier studies of 
dark energy show that these cosmic parameters obey power-law form of the time parameter. Next we plotted these 
parameters for various values of the BD parameter in different cosmic epochs. 

In figures 1 and 2, the scale factor is plotted is shown for different times. The positive and negative subscripts of 
a{t) correspond to two values given in (29). The behavior of scalar field (j) is plotted in figures 3 and 4. It is shown 
that the positive component of the field 0+ increases faster for small values of uj and vice versa for (/)_ . In figures 5 
and 6, the expressions for f+ and /_ are plotted. It is observed that these functions behave more like exponential 
functions which shows that the coupling between the BD field and matter increases with time. In figures 7 and 8, we 
provide the behavior of the BD potential which is increasing (decreasing) for V+ {V-) against time for large (small) 
values of the BD parameter. Figure 9 gives the variation of the BD potential against the BD field while the figure 10 
provides the variation of the scale factor in the limit when the BD parameters vanishes (becomes negligible). 

In figures 11, 12 and 13, we have plotted statefinder parameters and the deceleration parameter of our model against 
the BD parameter. It is observed that deceleration parameter is positive, thereby producing a decelerated Universe. 
Though the model suggests the expansion of the Universe, yet the model does not explain cosmic acceleration. The 
stability analysis of the solutions is important. But it will be beyond the scope of our paper. Such analysis can be 
done using the similar methodology, as it has been done in (TB], by defining dimensionless variables. The stability of 
BD with Chameleon field has been done already in [l7| for the same system of equations as we have. 
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